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Abstract 

We report on the experimental investigation of gluing bifurcations in the analog electronic circuit which models a 
fT^ Idynamical system of the third order: Lorenz equations with an additional quadratic nonlinearity. Variation of one of 
T-H the resistances in the circuit changes the coefficient at this nonlinearity and enables transition from the Lorenz route to 
^^ chaos to a different scenario which leads, through the sequence of homoclinic bifurcations, from periodic oscillations of 



Oh 



the voltage to the chaotic state. A single bifurcation "glues" in the phase space two stable periodic orbits and creates a 
'new one, with the doubled length: a bifurcation sequence results in the birth of the chaotic attractor. 
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.1. Introduction 

' Studies of routes from order to chaos in families of 
low-dimensional dynamical systems usually started with 
jtheoretical and numerical investigations, which were soon 
followed by detailed experiments. Accurate experimental 
verification not only recovered the qualitative picture of 
•the transition but was also able to resolve its quantitative 
characteristics. Universal sequence of period-doubling bi- 
furcations, onset of chaos via the breakup of the quasiperi- 
pdic oscillations as well as scaling laws in different kinds of 
intermittency have been documented in numerous mechan- 
ical, hydrodynamical, optical, chemical etc. experiments 
■(see references to original publications e.g. in [l|, l3|). 

There is, however, a group of scenarios which is well 
understood theoretically, but has received less attention 
from experimentalists: these are the routes to chaos via 
the sequences of so-called gluing bifurcations. Along such 
routes, pairs of stable periodic orbits come close in the 
phase space, recombine and form new stable periodic or- 
bits which are more complicated than the original ones. 
Each recombination is mediated by two trajectories, ho- 
moclinic to the same saddle point. Homoclinic bifurca- 
tions often occur in nonlinear dynamics: in the Lorenz 
equations [3|, the birth of the chaotic set in the phase 
space owes to the so-called "homoclinic explosion" [J, |5[ . 
The gluing bifurcation differs from the homoclinic explo- 
sion both in the number of newborn periodic orbits and 
in their stability. The explosion generates, in a single act 
of creation, the countable set of periodic orbits, each of 
them asymptotically unstable; taken together, they form 



the kind of "skeleton" for the emerging chaotic attractor. 
In contrast, a gluing bifurcation, taken alone, produces 
just one or two stable periodic orbits 6]. However, in the 
course of the sequence ("scenario") of such bifurcations 
the shape of the attracting orbit gets more and more in- 
volved, and its length grows, until the whole development 
culminates in the chaotic attractor. Which of two cases - a 
homoclinic explosion or a gluing bifurcation - takes place 
in the particular family of dynamical systems, is entirely 
determined by the ratio of the two leading eigenvalues of 
the Jacobian computed at the saddle point. 

From the point of view of an experimentalist, gluing 
scenarios have an unpleasant feature: vulnerability of their 
basic building blocks. Unlike periodic orbits, saddle con- 
nections in generic dynamical systems are structurally un- 
stable. Furthermore, a gluing bifurcation assumes coexis- 
tence, at the same set of parameter values, of two homo- 
clinic trajectories to the same saddle point. Since, gener- 
ically, formation of the homoclinic trajectory to a saddle 
point of a dissipative system is a codimension-one event, 
every gluing bifurcation is a codimension-two phenomenon. 
Hence, the accomplishment of the gluing scenario requires 
either the perfect mirror (Z2) symmetry of the system 
(which transforms a homoclinic orbit into another one, so 
that they exist in pairs), or the ability to track in the pa- 
rameter space the sequences of codimension-two events. 
This makes gluing bifurcations a difficult object for labo- 
ratory studies: they are sensitive both to fluctuations and 
to imperfections of the experimental setup. Most occur- 
rences of gluing scenarios in families of dynamical systems 
were reported in theoretical and numerical studies: in the 
context of hydrodynamics (3, iSj iBjj nematic liquid crys- 
tals [lOl lll| and optothermal devices 12| . In the experi- 
ments, separate gluing bifurcations were identified in the 
Taylor-Couettc flow of the viscous fluid between two cylin- 
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ders [13[ and in the Chua oscillators [IJ]. Here, we report 
on our experimental investigation of the sequence of gluing 
bifurcations in an analog electronic circuit. 

2. Theoretical predictions 

2.1. General considerations 

In the phase space of a continuous dissipative dynam- 
ical system, the basic ingredient of the gluing bifurcation 
is an equilibrium of the saddle type with one-dimensional 
unstable manifold. We start from the situation when the 
system, like the Lorenz equations, possesses a symmetry 
which transforms into each other two components of this 
manifold. When the parameters of the system are var- 
ied, location of invariant manifolds in the phase space 
varies as well. At certain combinations of parameters, 
one of the components of the unstable manifold can re- 
turn back to the saddle along the stable manifold and 
form the homoclinic orbit. Symmetry ensures that the 
second component returns as well: homoclinic orbits come 
in pairs. Sequence of events (scenario) which accompany 
the birth/destruction of homoclinic orbits, depends on the 
leading eigenvalues of linearization of the equations at the 
saddle point. Since unstable manifold is one-dimensional, 
there is just one positive eigenvalue, denoted below as A+. 
Below we restrict ourselves to the case when the closest 
to zero negative eigenvalue A_ is real. The "saddle in- 
dex" v — I A_ I A+ indicates which of the two properties - 
contraction or expansion - dominates the phase space in 
the neighborhood of the fixed point, and, thereby, gov- 
erns the stability of the bifurcating solutions. In absence 
of symmetry, destruction of a single homoclinic trajectory 
creates a unique periodic orbit which is stable if i^ > 1 and 
unstable otherwise [15[. Presence of the second, symmet- 
ric homoclinic trajectory enriches dynamics: all trajecto- 
ries which leave the vicinity of the saddle, are re-injected 
back. In this case, under v < 1 the countable set of un- 
stable periodic orbits, as well as a continuum of recurrent 
trajectories is simultaneously born from the pair of homo- 
clinic orbits [5|; this "homoclinic explosion" is a crucial 
step in the subsequent formation of the Lorenz attractor. 
The situation for z/ > 1 is simpler: here two stable peri- 
odic orbits approach the saddle point, and are "glued to- 
gether" forming two homoclinic orbits. When the pair of 
homoclinic trajectories breaks up, the new stable symmet- 
ric periodic orbit is left in the phase space: it is obtained by 
concatenation of the previously existing ones. The length 
and the number of loops (turns) of the attracting trajec- 
tory in the phase space is doubled, like in the case of the 
period-doubling bifurcation. Notably, in contrast to the 
period doubling bifurcation, the temporal period displays 
unbounded growth when the system approaches the bifur- 
cation: the period becomes infinite when the homoclinic 
trajectories are formed. 

Further variation of parameter can result in the se- 
quence of secondary gluing bifurcations: unstable mani- 
fold can return to the saddle after performing several turns 



in the phase space. Before this, a symmetry-breaking bi- 
furcation should take place: the newborn symmetric pe- 
riodic orbit is destabilized in the course of the pitchfork 
bifurcation, and two mutually symmetric orbits branch 
from it. These two orbits approach the unstable mani- 
fold of the saddle and coalesce in the next gluing bifur- 
cation. As a consequence, the new stable periodic orbit 
is born, which has four times more loops than the origi- 
nal ones. The subsequent scenario consists of alternating 
gluing- and symmetry-breaking bifurcations that eventu- 
ally end in the formation of the chaotic attractor which 
has a two-lobe shape, reminiscent of the Lorenz attrac- 
tor. The sequence of the bifurcational values of the pa- 
rameter has been shown to converge at the exponential 
rate [IGIlZI- In contrast to the period-doubling scenario, 
this rate is not the unique universal constant: renormal- 
ization group analysis shows that the universality class is 
completely predetermined by the (in general, non-integer) 
value of the saddle index v > 1. 

Remarkably, the attractor which is formed in the course 
of this bifurcation scenario, occupies a certain intermediate 
position between order and chaos: the Fourier spectrum of 
the trajectory is neither discrete like in case of regular dy- 
namics, nor absolutely continuous like in case of a chaotic 
or stochastic process, but is supported by the fractal set. 
Accordingly, observables are characterized through long- 
range non-exponential correlations [13] ■ 

If the symmetry between the components of the un- 
stable manifold is violated, each gluing bifurcation is a 
codimension-two event. On the plane of two parameters, 
there are numerous paths which lead from order to chaos 
via the formations of secondary homoclinic orbits; each of 
these paths is characterized by its own scaling constants [18J, 

iia,i2fl[2i|. 



2.2. The model set of equations 

While choosing a dynamical system which should serve 
as a candidate for experimental modeling by an electronic 
circuit, it is reasonable to minimize both the order of 
the system and the complexity of its linear and nonlin- 
ear terms. A unique gluing bifurcation can take place 
on the phase plane, but for occurrence of sequences of 
such bifurcations the phase space should be at least three- 
dimensional. A natural candidate is a system of the Lorenz 
equations [3]: they possess the desired symmetry and are 
capable to exhibit homoclinic bifurcations. In the three- 
dimensional parameter space of the Lorenz equations, for- 
mation of two symmetric homoclinic orbits has codimen- 
sion one: it takes place upon the two-dimensional surface. 
In the phase space, a homoclinic bifurcation is a nonlo- 
cal event. For this reason, the values of parameters under 
which the homoclinic orbits in dissipative dynamical sys- 
tems are formed, typically cannot be explicitly obtained in 
the closed form and should be found numerically. There 
are several analytical techniques which allow to exclude 
presence of homoclinic orbits and/or provide bounds for 
their existence in the parameter space 22], |23|, |2J, |25[. 



We are, however, not aware of rigorous theoretical resuhs 
which would allow to conclude on the location of the whole 
bifurcation surface in the parameter space of the Lorenz 
equations. Numerical evidence indicates that the two- 
dimensional surface of homoclinic bifurcation lies in the 
part of the parameter space which corresponds to i' < 1; 
hence, the "canonical" Lorenz equations display no gluing 
bifurcations. To circumvent this obstacle, we introduce an 
additional term and consider the set of equations 



X = a{y — x) ^^ Ay z 
y = Rx — y — X z 
z = xy — b z 



(1) 



Here, cr, R and b are the conventional Lorenz parameters, 
whereas A parameterizes the added nonlinear term in the 
first equation. At ^=0 Eqs (P) turn into the Lorenz equa- 
tions. The system ^ is reminiscent of the model, em- 
ployed in [Tj] for studies of thermal convection in the layer 
of fluid subjected to high-frequency modulation of gravity. 
We fix the traditional values of the parameters a—\0 and 
6=8/3 [3|, and vary the remaining parameters R and A. 

Among the properties of equations ([T]) we list only 
those which are relevant for gluing bifurcations and for 
experimental modeling. Equations ([Ij are invariant with 
respect to the transformation {x — >■ — x, y — >■ — y}. The 
origin x = y = z = Q is the equilibrium which is stable 
for R < 1 and is a saddle with one-dimensional unstable 
manifold in the parameter range R > 1. At R — \ the 
pitchfork bifurcation takes place at the origin; this bifur- 
cation is supercritical for A < a and subcritical otherwiso 
Like in the original Lorenz equations, variation of parame- 
ters can produce pairs of structurally unstable homoclinic 
orbits to the saddle point; these orbits leave the vicinity 
of the origin along the xy-plane and return to it along 
the z-axis; under a = 10, b = 8/3, A = Q this happens at 
R = 13.926 The saddle index of the origin equals 
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-CT - 1 + ^J{a -ly + AaR 



(2) 



Since the value of v is ^-independent, variation of A al- 
lows to study the effects of additional nonlinearity under 
constant saddle index. The value of v is smaller than 1 
for R > Ri^ = {b + <t)(6 + l)/cr and exceeds 1 otherwise. 
Numerical integration (we have used the recurrent Taylor 
algorithm of the 30th order with variable time-step) in- 
dicates that increase of A lowers the critical value i?hom, 
required for the formation of the pair of homoclinic or- 
bits. Therefore, at small positive values of A breakup of 
homoclinic orbits is followed by the homoclinic explosion 
and the Lorenz scenario of transition to chaos, whereas the 
sufficiently large values of A ensure the inequality v > 1, 
so that the sequence of gluing bifurcations is observed. 



'^In the latter case, the pitchfork is preceded by the saddle-node 
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3. Experimental setup 

For our measurements we took the electronic circuit 
with the help of which the authors of [26| reproduced the 
dynamics on the Lorenz attractor. In order to take ac- 
count of the additional nonlinear term in the first equation 
of (H)), we introduced an additional analog multiplier (the 
rightmost multiplier in Fig. [1]). 




Figure 1: Basic electronic circuit which models Eqs |^. Nominal 
values of resistors and capacitors are shown on the scheme. The 
analog multipliers (denoted by crosses) are of the type AD633AN. 
Variables «, i), w are output voltage values of operational amplifiers 
UIB, U2B and USB, respectively 

The cited values of electrical characteristics correspond 
to fixed values tr = 10 and b = 8/3 in Eqs (P). Changing of 
parameters A and R is enabled by variation of respective 
resistances R17 and R5. The value of A equals 10 divided 
by the value of R17 in fcfi; the value of R equals 0.01 of 
the resistance R5 in \Xl. 

To ensure that all voltages in the circuit stay inside the 
operating range of dynamical multipliers (from -10 to 10 
V) , the original dependent and independent variables need 
to be rescaled: we assume that the dimensionless potential 
differences a;, y and z are measured in volts, and proceed 
to u = cc/S, V = y/5, w = z/\0. For the time measured 
in "seconds", we introduce r = t/T with r=100. This 
recasts Eqs ([T]) into 



= T{a{v-u) + lQAvw) 
= T (Ru — V — 10 uw) 
~ T (2.5 uv — bw) 



(3) 



bifurcation which occurs at R - 



where prime denotes differentiation with respect to r. Fi- 
nally, it should be noted that, due to inevitable presence 
of nonlinearities in the multipliers, the real circuit, in fact, 
does not reproduce the perfect symmetry of Eq.Q with 



respect to the simultaneous change of sign of voltages u 
and V. 

In order to initiate the trajectories close to the unstable 
manifold of the saddle point, the circuit was grounded for a 
short time. The subsequent growth of voltage began from 
rather small values. 

Voltages were measured and recorded at a rate of 5000 
records per second. We did not apply any kind of filtering; 
stability of the scheme and smallness of time-step allowed 
us to plot quite smooth phase portraits of the system di- 
rectly from the experimental data. 

4. Results 

4-.1. Lorenz scenario of transition to chaos 

We measured and recorded time-dependent voltages in 
the circuit for various values of the parameters A and R 
at fixed cr=10 and 6=8/3. The equilibrium u—v=w=Q is 
stable for R <1 and is a saddle-point in the range R>1. 
Measurements on the circuit and numerical integration of 
Eqs ll]) confirm existence of the curve i?hom(^) upon which 
principal homoclinic orbits to this saddle are formed. Both 
in the experiment and in numerical studies this event is 
marked by change in the asymptotics of large t for trajec- 
tories with initial conditions near the origin. Prior to the 
bifurcation, the trajectory which starts with small positive 
values of u and v tends to the attractor (equilibrium or os- 
cillatory state) which is contained in the half-space u > 0. 
After the bifurcation, this trajectory makes a loop in the 
half-space u > 0, crosses into the opposite half-space w < 
and settles on the attractor there. 

Under low values of A, formation of homoclinic orbit 
occurs at t^ < 1; this corresponds to the Lorenz scenario 
of transition to chaos. Immediately beyond this curve 
the chaotic set, created by the homoclinic explosion at 
-Rhom(^), is unstable; its presence is reflected in the ex- 
istence of chaotic transients which precede relaxation to 
stable equilibria ( "metastable chaos" [221 )■ On the param- 
eter plane, the domain corresponding to this kind of dy- 
namics is bounded from above by the curve i?ch {A) which 
marks stabilization of the chaotic set and emergence of the 
chaotic attractor. Upon this curve the unstable manifold 
of the saddle point hits the stable manifold of the unsta- 
ble periodic orbit i5||. Numerically, we detected the latter 
curve by tracking the separatrices of the saddle point; in 
the experiments Rch{A) was identified as the lowest value 
of R (or the leftmost value of A) where the circuit pos- 
sessed non-decaying chaotic dynamics. For values of R be- 
yond i?ch(^), we observed irregular oscillations of voltages, 
with phase portraits strongly reminiscent of the Lorenz at- 
tractor (Fig.©. 

4-2. Sequences of gluing bifurcations 

According to the above reasoning, homoclinic explosion 
should be replaced by the gluing bifurcation when the cor- 
responding value of i?hom gets below i?^ =209/45=4.6444. . .; 
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Figure 2: Lorenz-like dynamics aX v <\ with cr=10 and 6=8/3. 
(a); Bifurcations on the parameter plane. SoUd line iJhom (numer- 
ics), and crosses (experiment): homocUnic explosion. Dashed line 
-Rch (numerics) and circles (experiment): onset of chaotic motion, 
(b): Projection of the phase portrait at ^4=1.5, -R=3.6 (experiment). 



this corresponds to the range A > 4. Indeed, we observed 
in the experiment the Lorenz-like chaotic attractors for 
A < 4 and gluing bifurcations for A > 4. 

We illustrate the transformation of the phase portrait 
along the gluing bifurcation scenario with the plots of pro- 
jections which correspond to increase of A at constant val- 
ues of other parameters. 
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Figure 3: Principal gluing bifurcation in the circuit: Evolution of the 
attracting set. Parameter values: (t=10, 6=8/3, iJ=3.6. 
(a): j4=6.4; two stable periodic orbits, (b): ^4=6.55; two stable 
periodic orbits, (c): A=I.\ formation of two trajectories, homoclinic 
to the equilibrium, (d): A=7.05; stable self-symmetric orbit with 
two loops. 

By varying initial conditions, we are able to identify in 
the phase space of the circuit two attracting closed trajec- 
tories which are almost symmetric to each other (Figj^Ji'). 
As the parameter A is increased, these two orbits come 
closer (FiglSb), approach the invariant manifolds of the 
saddle equilibrium at the origin and form the pair of sym- 
metric homoclinic orbits to this saddle (Fig 13]:). A slight 
increase of the parameter A makes the homoclinic orbits 
break up and disappear; the only attractor of the sys- 



tern, plotted in Fig|3jl, is the periodic orbit with two loops 
which is (up to instrumental resolution) invariant under 
the symmetry transformation {u — )■ —u, v — >■ —v} . Close 
to the homoclinic bifurcation, the measured values of the 
period Tq of the oscillations match the known logarithmic 
asymptotics: To{A) - - log(Ahom - A). 

As seen in Fig|4j the oscillations are strongly anhar- 
monic: the overwhelming part of the period is spent in 
nearly motionless state. This corresponds to long hover- 
ing of the trajectory in the vicinity of the saddle point. 
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Figure 4: Temporal evolution close to the gluing bifurcation, 
(a): A=7.0; oscillations on limit cycles with 1 loop prior to the bi- 
furcation, (b): A=7.05; oscillations on the limit cycle with 2 loops. 
Other parameters like in Fig. \3\ 

Further stages of the evolution of the attracting trajec- 
tory are sketched in FigO The symmetric limit cycle with 
two turns (Fig[3Jl) loses stability as a result of the pitch- 
fork bifurcation. One of the two resulting stable asym- 
metric limit cycles is shown in Fig|5^. As the parameter 
is further increased, these two limit cycles approach the 
invariant manifolds of the equilibrium, and the secondary 
gluing bifurcation takes place: two homoclinic orbits with 
2 turns are formed. Their subsequent breakup produces a 
single symmetric stable orbit with 4 turns (FiglSj)). Fur- 
ther, the events are repeated on the new level: the sym- 
metric orbit with 4 turns is destabilized, and two stable 
asymmetric ones are born (Fig|SJ;) . These asymmetric pe- 
riodic orbits merge in the next gluing bifurcation, form a 
pair of homoclinic orbits with 4 turns, and leave the stable 
periodic orbit with 8 turns, shown in FiglSJi. 

In the perfectly symmetric setup this bifurcation se- 
quence would continue, each time doubling the number of 
turns of the attracting periodic orbit in the phase space [16| . 
Presence of asymmetry is known to disable the complete 
sequence of homoclinic "doublings" 28|, |29j ; mathemati- 
cally, this owes to the fact, that linearization of the cor- 
responding renormalization operator near its fixed point 
possesses an additional unstable direction, responsible for 
the asymmetry. Finite resolution of our measurements (fi- 
nite step-size for variable resistance R17 in Fig.[T]) as well 
as the inevitable asymmetry in the electronic circuit did 
not allow us to resolve the further stages of the gluing pro- 
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Figure 5: Further stages of the gluing scenario (parameters: cr=10, 
b=8/3, -R=3.6). (a): A=8.05; asymmetric stable periodic orbit with 
two turns (b): j4=8.10; symmetric stable periodic orbit with 4 turns, 
born from the secondary gluing bifurcation, (c): 71=9.37; asymmet- 
ric stable periodic orbit with 4 turns; (d): ^4=9. 38; symmetric stable 
periodic orbit with 8 turns, born from the third gluing bifurcation. 



cess: instead we observe rapid emergence of the chaotic 
attractor (top row of Fig. [6]). 

A characteristic feature of chaotic oscillations in this 
circuit, well recognizable in the oscillograms of Fig.[Sl^a,e), 
is presence of relatively long nearly quiescent plateaus: 
they correspond to hovering near the saddle point. Fur- 
ther increase of the parameter A has disclosed that the 
range in which oscillations are chaotic, is interspersed by 
numerous windows in which various stable symmetric and 
asymmetric periodic oscillations are observed. One of such 
periodic states is shown in Fig. [nijc,d). In terms of the 
parameter A, periodic windows are bounded from above 
by secondary homoclinic bifurcations and serve as start- 
ing states of further gluing scenarios. An example of bi- 
furcation, in which two periodic regimes with three turns 
glue together and create a periodic orbit with six turns 
is shown in Fig. Etgjh)- It should be noted, that increase 
of the resistance R17 enhances asymmetry in the circuit, 
therefore, in general, homoclinic trajectories formed, re- 
spectively, by the "left" and the "right" components of 
the unstable manifold occur at slightly different values of 
the parameter A. 

5. Discussion and Outlook 

We have demonstrated in the experiment on the analog 
electronic circuit the Lorenz scenario of transition to chaos 
as well as several initial stages of the sequence of gluing 
bifurcations. The circuit mimics the modified Lorenz sys- 




Figure 6: Experimental records of dynamics beyond the first gluing 
sequence. (a,b): Chaotic attractor at A=9.39. (c,d): Limit cycle 
with 5 turns at A=10.17. (e,f): Chaotic attractor at j4=13.46. (g,h) 
gluing bifurcation, mediated by homoclinic orbits with 3 turns at 
A=11.51. Other parameters: a=lO, b=8/3, R=3.6. 



tern, and we utilized the fact that modification allows, by 
variation of a single resistance, to move across the param- 
eter space without changing the value of the saddle index. 
Observed sequence of transformations of phase portraits 
provides an unambiguous qualitative confirmation of the- 
oretical predictions. For a quantitative comparison (con- 
vergence rate of the bifurcation scenario and other scaling 
constants), the experimental resolution - in particular the 
step-size of the active parameter - should be improved. 

In a wider context of dynamical systems with gluing bi- 
furcations, it should be noted that a straightforward cod- 
ing procedure maps temporal evolution of voltages onto 
the binary alphabet: a symbol 1 is assigned to each turn 
of the orbit in the half-space m > 0, and a symbol marks 
each turn in the half-space u < 0. In this sense, a pair 
of binary codes which correspond to two components of 
the unstable manifold of the saddle, provides a complete 
characterization of dynamics [19| . In presence of symme- 



try, one code is sufficient; remarkably, it is related to the 
kneading sequence of the family of logistic mappings. In 
case of noticeable asymmetry, both kneading sequences are 
needed; here, certain aspects of the description are analo- 
gous to the formalism for (in general, discontinuous) cir- 
cle mappings: rotation numbers, devil staircase etc. |21| . 
Since our experimental results match well the theoretical 
description of the route to chaos for symmetric systems, we 
may hope that explicit introduction of controlled asymme- 
try into equations ^ and into the analog circuit of Fig. [T] 
should allow to detect in the circuit bifurcation also those 
scenarios which theory predicts for asymmetric gluing bi- 
furcations. 
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